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I.  INTRODUCTION 


To  accurately  characterize  dynamic  models  of  complex  structures, 
experimental  measurement  of  the  structure's  natural  modes  of  vibration  is 
required.  Typically,  the  measurement  of  these  modes  is  a  costly  and 
technically  demanding  taste  that  requires  both  experimental  and  analytical 
expertise.  Over  the  years,  methodologies  have  been  developed  to  improve  the 
accuracy  of  measured  modes  and  reduce  the  complexity,  and  therefore  the  cost, 
of  the  experimental  procedures.  Generally,  these  methodologies  can  be  cate¬ 
gorized  into  two  groups.  The  first  group  consists  of  procedures  that  attempt 
to  establish  natural  modes  of  vibration  by  direct  measurement  of  the  test 
article's  forced  vibration.  The  second  group  consists  of  those  procedures 
that  attempt  to  identify  natural  modes  of  vibration  by  post-test  analysis  of 
frequency  response  functions  or  post-test  analysis  of  time  domain  free  re¬ 
sponse  data. 

To  date,  procedures  that  attempt  to  establish  natural  modes  of  vibration 
by  direct  measurement  use  one  or  more  shakers  to  exert  sinusoidally  varying 
forces  on  the  test  article.  The  frequency  of  excitation  and  the  relative 
force  levels  of  the  multiple  shakers  are  adjusted  to  best  isolate,  from  all 
other  modes,  the  target  mode  response.  The  mode  parameters  are  then  estab¬ 
lished  from  direct  measurement  of  the  forced  vibration. 

Over  the  years,  test  methodologies  for  adjusting  the  relative  force 
levels  of  the  multiple  shakers  have  been  proposed  (e.g.,  Refs.  1-3).  These 


^Lewis,  R.  C.  and  Wrisley,  D.  L. ,  "A  System  for  the  Excitation  of  Pure 
Natural  Modes  of  Complex  Structures,"  Journal  of  Aeronautical  Sciences, 

Nov.  1950,  pp.  705-722. 

^Aaher,  G.  W. ,  "A  Method  of  Normal  Mode  Excitation  Utilizing  Admittance 
Measurements,"  IAS  Proceedings  of  the  National  Specialists  Meeting — Dynamics 
and  Aeroelasticity ,  Fort  Worth,  Texas,  6-7  Nov.  1958,  pp.  69-76. 

^Craig,  R.  R. ,  Jr.  and  Su,  Y.  W.  T. ,  "On  Multiple-Shaker  Resonance 
Testing,"  AIAA  Journal,  Vol.  12,  No.  7,  July  1974. 


procedures  have  relied  on  phase  coherence  of  response  as  Che  criterion  for 
successful  isolation  of  a  target  mode.  For  example,  the  methodology 

introduced  by  Lewis  and  Wrisley^  requires  that  the  shakers  oe  brought  into 
play  one  at  a  time.  As  each  new  shaker  is  introduced,  all  the  shaker  force 
levels  are  readjusted  in  an  attempt  to  obtain  a  force  distribution  in  phase 
with  response  velocity,  and  thus  establish  displacement  responses  at  all 
degrees  of  freedom  in  phase  with  each  ocher.  More  recently,  Anderson"* 

proposed  a  test  procedure  Chat  does  not  rely  on  phase  coherence  of  response  as 
a  criterion  for  mode  isolation.  Instead,  the  procedure  relies  on  maximizing  a 
target  mode  quadrature  response  relative  to  the  quadrature  responses  of  modes 
"close"  in  frequency  to  Che  target  mode.  This  is  accomplished  by  using 

successive  measured  estimates  of  shaker  location  mode  shape  values  to 
establish  a  target  mode  excitation  vector  nearly  orthogonal  to  modes  "close" 
in  frequency  to  the  target  mode. 

Several  investigators  have  proposed  test  and  analysis  procedures  chat 

attempt  to  establish  natural  modes  of  vibration  by  post-test  analysis  of 
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frequency  response  functions.  Morosow  and  Ayre  and  Stable  have  proposed 
analysis  procedures  to  extract  mode  shapes  from  measured  response  to 
nonselective  sinusoidal  excitation.  More  recently,  the  introduction  of  the 
Fast  Fourier  Transform  (FFT)  algorithm  and  the  availabililty  of  digital  data 
processing  equipment  have  led  to  Che  development  of  test  procedures  that 
make  no  attempt  to  cause  the  test  article  to  vibrate  at  discrete  frequencies. 
Instead,  all  modes  within  Che  frequency  range  of  interest  are  simultaneously 

^Kennedy,  C.  C.  and  Pancu,  C.  D.  P.,  "Use  of  Vectors  in  Vibration 
Measurement  and  Analysis,"  Journal  of  the  Aeronautical  Sciences,  Vol.  14, 

No.  11,  Nov.  1947. 

^Anderson,  J.  E.,  "Another  Look  at  Sine-Dwell  Mode  Testing,"  AIAA 
Journal  of  Guidance,  Control  and  Dynamics,  Vol.  5,  No.  4,  July-August  1982, 
pp.  358-365. 

^Morosow,  G.  and  Ayre,  R.  S.,  "Force  Apportioning  for  Modal  Vibration 
Testing  Using  Incomplete  Excitation,"  Shock  and  Vibration  Builecin  No.  48, 

PC.  1,  1978. 

^Stable,  C.  V.,  Jr.,  "Phase  Separation  Technique  for  Ground  Vibracion 
Testing,"  Aerospace  Engineering,  July  1962. 


excited  with  either  a  single  broad-band  randomly  varying  force,  multiple 
uncorrelated  broad-band  randomly  varying  forces,  or  an  impulsive  force.  The 
response  data  are  then  digitally  processed  into  estimates  of  steady-state 
frequency  response  functions,  from  which  modal  parameters  are  estimated  using 
any  one  of  several  candidate  procedures  (e.g.,  Refs.  8-10). 


In  addition  to  the  frequency  domain  identification  procedures,  Ibrahim 
and  Mikulick  *  have  proposed  analysis  procedures  that  attempt  to 
establish  natural  modes  of  vibration  by  post-test  analysis  of  measured  time 
domain  free  response  data.  Although  several  technical  papers  discussing  these 
procedures  have  been  published  (e.g.,  Refs.  11,  12,  and  13),  test  problems 
presented  to  date  have  emphasized  the  identification  of  natural  frequencies 
and  damping  values.  The  procedures'  ability  to  identify  mode  shapes  of  complex 
structures  still  needs  to  be  demonstrated. 


®Coppolino,  R.  N. ,  "A  Simultaneous  Frequency  Domain  Technique  for 
Estimation  of  Modal  Parameters  from  Measured  Data,"  SAE  Technical  Paper 
Series  811046,  Aerospace  Congress  and  Exposition,  Anaheim,  Ca.,  5-8  Oct.  1981. 

^Richardson,  M.  and  Potter,  R. ,  "Identification  of  the  Modal  Parameters 
of  an  Elastic  Structure  from  Measured  Transfer  Function  Data,"  20th 
International  Instrumentation  Symposium,  ISA  ASI  74250,  21-23  May  1974, 
pp.  239-246. 

lORichardson,  M. ,  "Modal  Analysis  Using  Digital  Test  Systems,"  Seminar 
on  Understanding  Digital  Control  and  Analysis  in  Vibration  Test  Systems,  A 
Publication  of  The  Shock  and  Vibration  Information  Center,  Naval  Research 
Laboratory,  Washington,  D.C.,  Part  2  of  2  Parts,  May  1975,  pp.  43-64. 

llIbrahim,  S.R.  and  E.C.  Mikulcik,  "A  Method  for  the  Direct  Identification 
of  Vibration  Parameters  from  the  Free  Response,"  Shock  and  Vibration 
Bulletin,  Bulletin  47,  Part  4,  Sept.  1977. 

Ibrahim,  S.R.,  "Modal  Confidence  Factor  in  Vibration  Testing,"  Journal 
of  Spacecraft  and  Rockets,  Vol.  15,  No.  5,  Sept.  1978. 

13Pappa,  R.S.,  "Close-Mode  Identification  Performance  of  the  ITD  Algorithm," 
AIAA/ASME/ASCE/AHS  24th  Structures,  Structural  Dynamics  and  Materials 
Conference,  Part  2,  2-4  May  1983. 


Ia  Che  past  few  years,  several  mode  survey  cescs  on  complex  structures 

have  beea  performed  such  that  selected  modes  were  measured  using  both 

mulci-shaker  sine-dwell  and  single-point-random  test  techniques.  With  mode 
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orthogonality  used  as  a  comparison  criterion,  published  evidence  ’  and 
ocher  unpublished  test  results  indicate  chat,  for  complex  structures, 
multi-shaker  sine-dwell  test  techniques  presently  yield  better  results  than 
single-point-random  test  and  mode  identification  techniques.  An  advantage 
presently  enjoyed  by  the  multi-shaker  sine-dwell  test  technique  is  that  the 

modes  of  vibration  are  established  before  disassembly  of  the  test  setup,  and 
any  deficiencies  in  the  quality  of  the  data  (as  caused  by  a  malfunctioning 
accelerometer,  for  example)  can  be  reconciled.  With  the  post-test  mode 
identification  procedures,  it  is  more  difficult  to  identify  "bad"  data, 
because  very  often  the  modes  of  vibration  are  not  available  until  months  after 
the  test,  and  no  additional  data  can  be  collected  because  schedule 
considerations  have  forced  the  disassembly  of  the  test  setup.  Another 

advantage  enjoyed  by  the  multi-shaker  sine-dwell  test  technique  is  that  the 

modal  parameters  are  the  result  of  direct  observations  of  the  test  article's 

forced  vibration.  With  the  random  or  impulsive  test  and  mode  identification 
techniques,  the  modal  parameters  are  established  from  "measured"  frequency 
response  functions  as  "best  estimates."  Therefore,  they  are  dependent  on  the 
numerical  techniques  (and  their  inherent  assumptions  and  limitations)  used  in 
the  identification  process. 

As  discussed  above,  mode  survey  test  procedures  can  be  catagorized 
according  to  whether  they  attempt  to  establish  natural  modes  of  vibration  by 
direct  measurement  or  attempt  to  identify  modes  of  vibration  by  post-test 
analysis  of  response  to  nonselective  excitation.  Procedures  that  establish 

l^Knauer,  C.  0.,  Peterson,  A.  J.,  and  Rendahl,  W.  8.,  "Space  Vehicle 

experimental  Modal  Definition  Using  Transfer  Function  Techniques,"  Presented 
at  the  SAE  National  Aerospace  Engineering  and  Manufacturing  Meeting,  Culver 
City/Los  Angeles,  Ca.,  17-20  Nov.  1975. 

^Perrante,  M. ,  Stable,  C.  V.  and  Breskman,  D.  G. ,  "Single-Point-Random  and 
Multi-Shaker  Sine  Spacecraft  Modal  Testing,"  Aerospace  Testing  Seminar, 

Los  Angeles,  Ca.,  Sept.  1979. 


natural  modes  of  vibration  by  direct  measurement  use  one  or  more  shakers  to 
exert  sinusoidally  varying  forces  on  the  test  article.  The  frequency  of 
excitation  and  the  relative  force  levels  of  the  shakers  are  adjusted  to  best 
isolate  the  target  mode  response.  Mode  parameters  are  then  established  from 
direct  measurement  of  the  forced  vibration. 

The  purpose  of  this  report  is  to  introduce  a  mode  survey  test  procedure 
that  uses  correlated  broad-band  randomly  varying  forces  to  accomplish  the  same 
task  that  until  now  required  harmonic  excitation  of  the  test  article.  The 
procedure  derives  from  the  recognition  that  response  data  of  a  structure, 
excited  at  multiple  locations  by  proportional  randomly  varying  forces,  ca  a 
analyzed  such  that  the  results  are  equivalent  to  those  obtained  had  he 
excitation  been  harmonic.  The  procedure  takes  advantage  of  multiple  s'  ar 
excitation  to  best  isolate  the  target  mode  response  from  all  other  m  .<? 
Mode  parameters  can  then  be  established  from  direct  observation  of  frequency 
response  functions. 


II.  FREQUENCY  RESPONSE  FUNCTIONS— A  REVIEW 


The  matrix  differencial  equation  of  motion  chat  descrioes  the  behavior 
of  a  structure  subjected  to  random  excitation  {A}f(t)  is  written  in  modal 
coordinates  as 

[  I  ]  {q  ( t  >  }  +  [256JJ  {q(t;  }  +  [u)2]{q(t)}  *  [<j>]  T{A>f  (t>  (1) 

where  the  coordinate  transformation  between  physical  coordinates  {X(tJ} 
and  modal  coordinates  (q(t)}  is  defined  by 

{X(t)}  -  [$Hq(t)}  ( 2 ) 

and  the  matrix  of  mode  shape  vectors  [<J>]  has  been  normalized  with  respect  to 
tne  system  mass  matrix  [M]  such  that 

[4>]T[M][4>1  -  [I]  (3) 


For  clarity  of  presentation,  we  shall  use  the  continuous  infinite-range 
Fourier  transform  to  establish  frequency  response  functions ,  Note,  however, 
that  we  shall  assume  later  that,  for  our  purposes,  frequency  response 
functions  obtained  with  digital  FFT  technology  are  equivalent  to  frequency 
response  functions  obtained  with  the  continuous  infinite-range  Fourier 
transform.  Proceeding,  we  take  the  Fourier  transform  of  Eq.  (1)  and  solve  for 
the  Fourier  transform  of  the  modal  coordinate  acceleration  vector  to  obtain 


{q(cu)  }  -  -u)2|[u)2-w2]2  +[2Cwnw]2yi^[u2-w2]-i[2Ca)n(u]J  [ 4>  1  T{a> f  C o>)  (4) 


ch 


where  Che  k  coordinate  equation  has  the  following  form: 


H  * 

£(u>)  - 


^  ♦ikAi 


Jt-1 


(8) 


factor  used  to  scale  random  excitation  time  history  at 
degree  of  freedom  l 

Fourier  transform  of  random  excitation  time  history  f(t) 


Taking  the  Fourier  transform  of  the  second  time  derivative  of  Eq.  (2)  and 
using  Eq.  (6)  to  substitute  for  Vt u>) ,  we  obtain  for  a  typical  physical 
coordinate  acceleration 

r 

xt«»>/fu> -2Xk*ian)?k  (») 

k«l 


The  right-hand  side  of  Eq.  (9)  is  recognised  to  be  the  response  function 
of  the  degree  of  freedom  of  a  multi-degree-of-f reedom  structure  that 
is  being  driven  at  multiple  locations  by  a  sinusoidally  varying  forcing 
function  of  frequency  «).  The  real  part  corresponds  to  the  coincident 
response  (i.e.,  the  component  of  acceleration  response  that  is  colinear  with 
the  reference  force  time  history),  and  the  imaginary  part  corresponds  to  the 
quadrature  response  (i.e.,  the  component  of  acceleration  response  that  is  90 
degrees  out  of  phase  with  the  reference  force  time  history).  Therefore,  Eq. 
(9)  implies  that  the  response  data  of  a  structure,  excited  at  multiple 
locations  by  proportional  randomly  varying  forces,  can  be  analyzed  such  that 
the  results  are  equivalent  to  those  obtained  if  the  structure  were  to  be 
excited  at  chose  same  locations  by  sinusoidally  varying  forces. 


III.  MODE  ISO'  \TION  LOGIC 


Assume  for  Che  moment  that  the  test  article  can  be  simultaneously  driven 
at  all  its  degrees  of  freedom  by  proportional  randomly  varying  forces. 
Equation  (9)  then  indicates  that  the  test  article  can  be  made  to  vibrate  in  a 
single  normal  mode  by  appropriately  adjusting  the  relative  shaker  force  levels 
such  that  Che  generalized  forces  are  zero  for  all  modes  except  the 
target  mode.  In  practice,  however,  Che  number  of  degrees  of  freedom  needed  to 
represent  a  complex  structure  greatly  exceeds  the  number  of  available 
shakers.  Therefore,  perfect  isolation  is  not  possible. 

Fortunately,  for  practical  purposes,  some  small  amount  of  contamination 
in  measured  modes  is  tolerable.  It  has  become  a  widely  accepted  practice  to 
judge  the  quality  of  measured  modes  by  their  mutual  orthogonality.  This  is 
accomplished  by  calculating  the  unit  normalized  generalized  mass  matrix  and 
comparing  the  magnitudes  of  the  off-diagonal  terms  to  a  predetermined  value. 
It  is  generally  accepted  in  the  industry  that  modes  of  acceptable  quality  have 
been  measured  if  these  off-diagonal  terms  are  no  greater  chan  0.10.  To  a 
first  approximation,  this  implies  that  a  target  mode  is  sufficiently  isolated 
for  measurement  if  the  amplitude*  of  Che  contaminating  mode  is  no  greater  than 
0.10  times  the  target  mode  amplitude,  i.e. 

<  (0.10)  Fc  (10) 

Therefore,  we  shall  proceed  with  Che  understanding  that  perfect  isolation  of  a 
target  mode  is  not  necessary  and  that  Che  number  of  shakers  need  not  equal  the 
number  of  degrees  of  freedom  in  the  test  article. 

It  was  noted  earlier  that  the  response  data  of  a  structure,  excited  at 
mulciple  locations  by  proportional  randomly  varying  forces,  can  be  analyzed 
such  that  Che  results  are  equivalent  to  chose  obtained  if  the  structure  were 


* 


Mode  shape  amplitude  refers  to  the  quadrature  component  of  the  total  response. 


Co  be  exciced  by  mule  ip Le  sinusoidally  varying  forces.  Therefore,  it  should 
be  possible  Co  adopc  Che  basic  principles  of  multi-shaker  sine-dwell  mode 
survey  cescing  and  apply  chem  co  mulci-shaker  random  mode  survey  Cescing. 

For  Che  multi-shaker  random  mode  survey  CesC  procedure  proposed  herein, 
Che  Anderson  mode  isolaCion  logic  of  seleccive  orChogonal  excitation 
(SOREX)^  has  been  adopted.  The  SOREX  logic  Cakes  advantage  of  Che  natural 
selectivity  of  lighcly  damped  sCructures  Co  greatly  amplify  cheir  forced 
responses  ac  resonant  frequencies.  A  consequence  of  chis  observation  is  recog¬ 
nition  Chat  careful  aCCenCion  Co  Che  selectivity  of  Che  applied  forces  general¬ 
ly  is  only  necessary  for  modes  in  close  frequency  proximity  Co  Che  targeC  mode. 
This  can  besc  be  illustrated  by  rewriting  Eq.  (10)  as 

(l£/l£)(F  /F  )  <  0.10  Ul) 

The  above  equation  states  that  adequate  target  mode  isolation  from  an  off- 
resonance  mode  will  exist  if  the  product  of  the  quadrature  admittance  ratio  and 
generalized  force  ratio  is  less  than  0.10. 

Equation  (11)  is  che  essence  of  Che  SOREX  procedure.  Basically,  the  pro- 

II  It 

cedure  recognizes  that  for  lighcly  damped  structures  the  ratio  (fi^/Q^)  will 
decrease  rapidly  with  increasing  frequency  separation  between  the  two  modes. 
Thus,  careful  attention  to  the  applied  forces,  which  determine  the  ratio 
(F^/Fg),  is  necessary  only  for  modes  in  close  frequency  proximity  to  the  target 
mode.  Data  presented  in  Ref.  5  and  test  experience  with  complex  spacecraft 
structures  indicate  that  if  a  mode  is  removed  in  frequency  from  the  target  mode 
by  store  than  about  10  to  15  percent,  it  can  generally  be  ignored  during 
measurement  of  Che  target  oxide.  Therefore,  the  applied  forces  will  generally 
be  required  only  to  isolate  Che  target  mode  responses  from  modes  "close"  in 
frequency.  To  accomplish  this,  the  SOREX  test  logic  is  proposed. 

Jiasically,  the  procedure  consists  of  using  successive  estimates  of  mode 
shape  values  to  establish  force  vectors  whose  corresponding  generalized  force 
for  che  target  mode  is  large  in  relation  to  the  generalized  forces  of  oxides 
"close"  in  frequency.  This  is  accomplished  by  establishing  isolation  groups 
chat  consist  of  one  or  more  target  modes  and  all  modes  within  about  10  to 
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15  percent  in  frequency.  To  obtain  optimum  isolation,  the  number  of 
excitation  shakers  must  equal  the  number  of  modes  in  the  isolation  group.* 
Typically,  for  complex  structures,  each  isolation  group  will  consist  of  2  to 
12  modes.  Therefore,  a  relatively  small  number  of  shakers  will  be  required. 

The  initial  shaker  force  levels  are  established  from  frequency  response 
functions  obtained  with  single  shaker  excitation.  The  response  functions  are 
used  to  estimate  the  shaker  location  mode  shape  values,  which  are  then  used  to 
calculate  the  relative  force  levels  to  be  applied  at  those  locations: 


A  A  A 

^1*  ^2*  **‘^m^ 


[7]uJ$1$2...$o]Vl 


where 


vector  of  forces  for  excitation  of  mode  at  gj. 


estimate  of  kth  mode  shape  vector  (quadrature 
response),  elements  of  shaker  position 
coordinates  only 

matrix  whose  kc^  column  represents  ideal  generalized 
force  vector  for  isolating  kc**  mode  from  other  modes 
in  isolation  group — typically  an  identity  matrix 


Each  calculated  force  vector  is  then  applied  to  the  structure,  the  frequency 
of  excitation  adjusted  to  maximize  quadrature  response,  and  new  values  of 
{^}  measured.  These  new  mode  shape  values  are  Chen  used  in  Eq.  (12)  to 
calculate  a  sec  of  refined  shaker  force  levels.  If  the  number  of  shakers  is 
equal  to  Che  number  of  modes  in  an  isolation  group,  which  we  shall  assume, 


*IC  is  possible  to  use  the  SOREX  test  logic  when  the  number  of  modes  in  an 
isolating  group  exceeds  the  available  shakers  (see  Ref.  5);  however,  the  iso¬ 
lation  will  generally  not  be  as  complete. 


•*« 


convergence  of  Che  force  vectors  to  those  -required  for  adequate  mode  isolation 
typically  occurs  within  one  or  two  iterations.  Once  isolation  is  achieved, 
the  entire  mode  shape  is  measured  and  its  orthogonality  checked  against  the 
ocher  measured  mode  shapes. 

The  test  procedure  proposed  by  Anderson  was  formulated  for  use  with 
harmonic  excitation  of  the  test  article.  However,  Eq.  (9)  indicates  that  if 
we  excite  the  structure  with  multiple  randomly  varying  forces,  the  response 
data  can  be  analyzed  such  Chat  Che  results  are  equivalent  to  Chose  obtained 
had  the  excitation  been  harmonic.  Therefore,  it  is  possible  Co  use  the  SOREX 
logic  to  establish  the  relative  magnitudes  of  multiple  randomly  varying  forces 
such  that  a  target  mode  is  isolated  for  measurement  from  modes  "close"  in 
frequency.  We  now  have  all  the  ingredients  needed  to  propose  a  multi-shaker 
random  mode  survey  test  procedure. 


IV.  MULTI-SHAKER  RANDOM  MODE  SURVEY  TEST  PROCEDURE 


The  basic  principles  of  multi-shaker  sine-dwell  mode  survey  testing  can 
be  combined  with  multi-shaker  broad-broad  random  excitation  to  establish  the 
following  procedure: 

1.  Use  single-point-random  excitation  to  establish  frequency  response 
functions  at  all  potential  shaker  positions  and  any  coordinates  whose  response 
will  aid  in  identifying  all  modes  in  the  frequency  range  of  interest. 


2.  Determine  preliminary  natural  frequencies  and  predominant  motions 
of  each  mode. 


3.  Measure  all  modes  which  are  removed  in  frequency  from  all  other 
modes  (i.e.,  0.85  <  uc  <  l«15  by  applying  a  randomly  varying  force  at  a 

location  on  the  test  article  that  exhibits  strong  response  in  the  mode  to  be 
measured.  Then,  process  the  response  data  into  frequency  response  func¬ 
tions.  Steps  11  and  12  define  how  the  modal  parameters  are  to  be  extracted 
from  the  response  functions. 


4.  Establish  isolation  groups  of  modes  sufficiently  close  in  frequency 
to  warrant  inclusion  in  determining  multiple  excitation  force  levels  (i.e., 
0.85  0>t  <  <  1.15  wt). 

5.  For  each  isolation  group,  select  for  shaker  locations  those 
coordinates  with  strong  resonances  and  phase  reversals  between  modes. 


6.  For  each  of  the  selected  coordinates,  extract  from  the  imaginary 
part  of  the  frequency  response  functions  established  in  step  1  the  peak  values 
associated  with  each  of  the  modes  in  an  isolation  group  and  form  [$]. 


-  «.  V'  -L  W  V"  *.  \  \  \ 
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A  A  A 

7.  Using  Eq.  (12),  calculate  [F^,  F2’**‘fnJ* 

8.  Apply  random  excication  at  the  coordinates  selected  in  step  5  such 
Chat  the  relative  force  levels  of  the  shakers  are  defined  for  each  target  mode 
by  {Ft>. 

A 

9.  For  each  {F^},  select  the  force  signal  from  one  shaker  as  a 
reference  signal,  and  process  the  response  data  from  each  shaker  location  into 
frequency  response  functions. 


10.  Review  frequency  response  functions  from  step  9  Co  determine  if 
adequate  isolation  has  been  achieved.  If  isolation  has  not  occurred,  form  a 
new  [$]  with  updated  shaker  location  mode  shape  values  from  step  9  and 
repeat  steps  7  through  10. 


11.  Record  the  mode  shape  as  the  maximum  values  of  the  imaginary  part 
of  the  frequency  response  function  of  each  coordinate.  The  natural  frequency 
is  recorded  as  the  frequency  at  which  the  imaginary  part  is  maximum. 


12. 

damping: 


where 


Use  the  real  part  of  the  frequency  response  function  to  calculate 


(13) 


Cc  “  critical  damping  ratio 

uib  ■  frequencies,  in  the  vicinity  of  the  target  mode 

frequency,  at  which  Che  real  part  of  the  frequency  response 
function  is  maximum  or  minimum 

bV  *  target  mode  natural  frequency 
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Repeat  steps  5  through  12  for  each  isolation  group. 


14.  Check  orthogonality  of  all  measured  modes.  If  unacceptable 
contamination  exists,  expand  the  isolation  group  of  the  contaminated  mode  to 
include  Che  contaminating  modes,  establish  a  refined  force  vector,  and 
remeasure  Che  contaminated  mode.  Note  that  with  Che  measured  mode  set,  a  good 
approximation  to  the  generalized  forces  can  be  calculated.  These  forces  can 
be  used  as  an  aid  in  determining  which  mode,  in  a  pair  of  modes  exhibiting 
poor  orthogonality,  should  be  remeasured. 


V.  DEMONSTRATION  OF  PROCEDURE  BY  NUMERICAL  SIMULATION 


To  demonstrate  Che  mu  Lei- shaker  random  mode  survey  test  procedure,  an 
analytical  simulation  of  a  test  was  performed  and  che  results  compared  to  (a) 
those  obtained  using  Che  multi-shaker  sine-dwell  mode  survey  test  procedure 
presented  in  Ref.  5,  and  (b)  the  eigenvectors  obtained  by  solving  che  syscem 
undamped  eigenvalue  problem.  The  analytical  model  employed  in  che  simulation 
had  the  following  dynamic  properties: 


[<»] 


1.954 

-0.742 

-2.360 

-0.209 

0.127 

0.009 

0.318 

-0.060 

-0.020 

1.959 

-2.340 

-0.763 

2.156 

-0.933 

2.103 

-0.205 

0.124 

0.008 

0.110 

0.029 

-0.005 

1.064 

-0.066 

2.975 

0.174 

0.189 

-0.001 

2.211 

2.116 

-0.752 

1.179 

2.922 

0.072 

-0.232 

-0.113 

0.008 

34.462 

r 

0.01 

36.098 

0.01 

o>  -  • 

37.503 

> 

(cl  -  < 

0.01 

■ 

121.179 

0.02 

166.086 

0.02 

.  356.219 

0.03 

The  response  of  Che  model  to  multiple  sinusoidally  varying  forces  was 
obtained  by  establishing  the  closed  form  complex  solution  to  che  equations  of 
motion.  The  real  part  of  the  solution  is  Che  coincident  response  and  che 
imaginary  part  of  the  solution  is  the  quadrature  response.  Tne  response  to 
multiple  randomly  varying  forces  was  established  by  numerically  integrating 
Che  equations  of  motion.  The  computed  response  time  histories  and  che 
reference  force  time  history  were  then  digitally  processed  into  frequency 
response  functions  from  which  the  modal  parameters  were  directly  inferred 
(without  curve  fits). 


The  test  simulation  was  initiated  by  exciting  the  analytical  test 


article  at  coordinate  X^  with  a  randomly  varying  force  and  processing  the 


response  data  at  selected  coordinates  into  frequency  response  functions.  From 
tnese  frequency  response  functions,  shown  in  Fig.  1*  for  the  frequency  range 
20-50  Hz,  it  was  determined  that  three  modes  were  present  in  the  frequency 
range  of  interest,  0  to  50  Hz.  Furthermore,  it  was  concluded  that  all  three 
modes  were  sufficiently  close  in  frequency  to  warrant  inclusion  in  the  same 
isolation  group. 


For  the  multi-shaker  excitation  phase  of  the  simulation,  coordinates 

X, ,  X0,  and  X.  were  selected  as  shaker  locations.  Estimates  of  mode 

1  3  4 

shape  values  at  these  coordinates  were  obtained  from  the  frequency  response 

functions  established  in  the  single-shaker  random  excitation  phase  of  the 

simulation.  The  mode  shape  values  were  taken  as  the  peak  imaginary 

response  at  the  estimated  natural  frequency  of  each  mode.  If  an  obvious  peak 

did  not  exist,  as  was  the  case  for  the  third  mode  at  coordinate  X.  (see 

4 

Fig.  1),  the  value  at  the  estimated  frequency  of  the  mode  was  used.  The 
estimates  of  mode  shape  values  were  then  used  in  Eq.  (10)  to  establish  the 
relative  excitation  force  amplitudes  needed  to  improve  the  isolation  of  each 


target  mode. 


The  acceleration  response  of  coordinate  X^  at  eacn  iteration  is 
shown  in  Fig.  2.  As  the  figure  illustrates,  the  isolation  of  each  target  mode 
improves  as  the  relative  excitation  force  levels  become  more  refined.  This  is 
especially  evident  for  the  second  mode,  where  it  can  be  observed  that  for 


single  shaker  excitation  the  dominant  response  of  coordinate  x^  was  that  of 


mode  1.  The  first  iteration  with  three  shakers  succeeded  in  isolating  the 
second  mode  from  the  first  mode;  however,  adequate  isolation  from  both  the 
first  and  third  modes  required  one  more  iteration. 


*For  ease  of  comparison,  all  frequency  response  functions  have  oeen  normalized 
to  provide  unit  generalized  mass  for  each  target  mode. 
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The  calculated  suaker  force  levels,  the  resulting  generalized  forces, 
and  the  unit  normalized  generalized  mass  matrix  obtained  at  each  iteration  are 
presented  in  Table  1.  As  the  table  illustrates,  the  single  shaker  excitation 
resulted  in  a  relatively  large  generalized  force  for  the  first  mode,  which 
partially  explains  the  dominant  response  of  this  mode  observed  at  coordinates 
X^,  X£»  X^,  and  X^  (see  Fig.  1).  Subsequently,  the  application  of 
multiple  shakers,  as  Table  1  also  illustrates,  resulted  in  the  target  mode 
generalized  force  being  maximized  relative  to  i.he  generalized  forces  of  the 
other  modes  in  the  isolation  group.  In  addition,  the  generalized  mass 
matrices  indicate  that  the  test  success  criterion  was  satisfied  after  the 
first  iteration:  all  off-diagonal  terms  in  the  unit  normalized  generalized 
mass  matrix  are  no  greater  than  0.10. 

The  mode  shapes  obtained  at  each  iteration  are  compared  in  Table  2  to 
the  mode  shapes  obtained  using  the  multi-shaker  sine-dwell  mode  survey  test 
procedure  presented  in  Ref.  5,  and  to  the  eigenvectors  obtained  by  solving  the 
system  undamped  eigenvalue  problem.  For  the  analytical  model  employed  in  the 
simulation,  mode  shapes  obtained  with  the  multi-shaker  random  excitation 
procedure  are  indistinguishaole  from  mode  shapes  established  with  the 
multi-shaker  sine-dwell  excitation  procedure  advocated  in  Ref.  5  (see 
Table  2).  In  addition,  excellent  agreement  exists  between  the  exact  mode 
vectors  and  those  established  in  the  simulations. 


fable  1.  Summary  of  Numerical  SimulaCion  Iterations 


SINE -DWELL 


ITERATION 


1.948 

0.318 

2.143 

0.111 

0.176 

1.211 


1.950 

0.319 

2.171 

0.110 

0.173 

1.157 


-0.374  -0.827 

-0.003  -0.061 

-0.574  -0.858 

0.049  0.029 

0.218  0.189 

3.079  2.922 


2.398  -2.361 

0.106  -0.020 
-1.506  2.103 

0.058  -0.005 

0.137  -0.001 

1.395  0.074 


1.955 

0.318 

2.157 

0.111 

0.174 

1.179 


-0.741 

-0.060 

-0.932 

0.030 

0.190 

2.923 


-2.358 

-0.020 

2.105 

-0.005 

-0.001 

0.065 


RANDOM 


ITERATION 


1.947 

0.316 

2.143 

0.110 

0.175 

1.215 


1.951 

0.318 

2.174 

0.110 

0.172 

1.150 


-0.309 

-0.856 

0.005 

-0.061 

-0.512 

-0.829 

0.053 

0.029 

0.220 

0.189 

3.096 

2.922 

2.381 

-2.361 

0.117 

-0.020 

-1.482 

2.103 

0.062 

-0.005 

0.152 

-0.001 

1.451 

0.075 

1.954 

0.317 

2.156 

0.110 

0.173 

1.178 


-0.741 

-0.060 

-0.932 

0.029 

0.189 

2.923 


-2.358 

-0.020 

2.106 

-0.005 

-0.001 

0.062 


1.954 

0.318 

2.156 

0.110 

0.174 

1.178 


-0.742 

-0.060 

-0.933 

0.029 

0.189 

2.922 


-2.360 

-0.020 

2.103 

-0.005 

-0.001 

0.072 


VI.  SOME  PRACTICAL  CONSIDERATIONS 


The  actual  performance  of  a  mode  survey  test  is  more  difficult  than  the 
mathematical  simulations  lead  us  to  believe.  Test  experience  indicates  that 
interaction  between  the  test  article  and  the  excitation  shakers  introduces 
phase  distortions  between  the  shaker  applied  forces.  For  harmonic  excitation, 
this  problem  is  overcome  by  delaying  in  time  one  shaker  force  signal  relative 
to  another  and  thus  adjusting  for  any  phase  distortions.  This  task  is 
considerably  more  difficult  for  random  excitation,  because  the  phase 
distortions  not  only  will  vary  from  shaker  to  shaker  but  also  will  vary  as  a 
function  of  frequency  at  each  shaker.  At  least  conceptually,  this  problem  can 
be  minimized  with  a  closed-loop  control  system.  This  added  complexity  needs 
to  be  considered  when  contemplating  a  multi-shaker  random  mode  survey  test. 

As  discussed  previously,  target  mode  isolation  from  modes  that  are  not 
included  in  the  isolation  group  is  primarily  due  to  separation  in  mode 
frequencies.  However,  occasionally  the  generalized  force  of  a  mode  not 
included  in  the  isolation  group  is  sufficiently  large  such  that  10  to  15 
percent  frequency  separation  from  Che  target  mode  is  not  sufficient.  Two 
corrective  actioas  can  be  taken.  The  simpler  approach  is  to  include  the 
contaminating  mode  in  Che  isolation  group,  deploy  an  additional  shaker,  and 
remeasure  the  contaminated  mode.  The  second  approach  is  to  move  the  shaker 
that  is  primarily  responsible  for  Che  large  response  of  the  contaminating 
mode.  First,  the  measured  target  mode  and  measured  contaminating  mode  must  be 
normalized  to  unit  generalized  mass.  Next,  the  mode  shape  values  at  the 
shaker  locations  should  be  compared.  Usually,  Che  contaminating  mode  will 
have,  at  one  of  the  shaker  locations,  a  considerably  larger  response  chan  Che 
target  mode.  The  shaker  attached  at  chat  location  should  be  redeployed, 
preferably  to  a  location  where  the  target  mode  response  is  greater  than  that 
of  Che  contaminating  mode. 


Although  numerical  simulations  can  serve  as  valuable  research  tools,  the 
true  test  o f  a  mode  survey  test  procedure  occurs  in  the  laboratory.  In  early 
1982,  the  proposed  multi-shaker  random  mode  survey  test  procedure  was 
successfully  used  to  establish  the  modes  of  vibration  of  a  spacecraft 
structure.  During  the  test,  several  modes  were  isolated  for  measurement  using 
multi-shaker  (three)  correlated  broad-band  random  excitation.  Each  shaker  was 
commanded  from  the  same  signal  generator;  however,  no  closed-loop  control 
system  was  used.  Of  the  ten  modes  measured,  eight  satisfied  the  test 
orthogonality  success  criterion  of  all  off-diagonal  terms  in  the  unit 
normalized  generalized  mass  matrix  being  no  greater  than  0.10.  The  one 
off-diagonal  term  greater  than  0.10,  which  was  0.19,  occurred  between  a  pair 
of  modes  in  which  the  resonant  components  were  inaccessible  to  the  shakers. 


VII.  CONCLUSIONS 


A  node  survey  test  procedure  has  been  introduced  Chat  comoines  random 
excitation  and  digital  processing  of  response  data  with  the  inode  isolation 
features  of  multiple  shaker  excitation.  The  proposed  procedure  uses  correlated 
broad-band  randomly  varying  forces  to  accomplish  the  same  task  that  up  to  now 
required  harmonic  excitation  of  the  test  article.  The  procedure  derives  from 
the  recognition  that  response  data  of  a  structure,  excited  at  multiple  loca¬ 
tions  by  proportional  randomly  varying  forces,  can  be  analyzed  such  that  the 
results  are  equivalent  to  those  obtained  had  the  excitation  been  harmonic. 
The  procedure  takes  advantage  of  multiple  shaker  excitation  to  best  isolate 
the  target  mode  response  from  all  other  modes.  Mode  parameters  are  then  esta¬ 
blished  from  direct  observations  of  frequency  response  functions. 

The  multi-shaker  random  mode  survey  test  procedure  has  been  demonstra¬ 
ted  herein  by  analytical  simulation  of  the  test.  The  results  of  the  simulation 
are  compared,  with  excellent  agreement  obtained,  to  values  established  with  a 
multi-shaker  sine-dwell  test  procedure  and  to  the  exact  results. 
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NOMENCLATURE 


veccor  of  elements  A^ 

factor  used  to  scale  random  excitation  time  history  at 
Ath  degree  of  freedom 

generalized  force 

vector  of  forces  for  excitation  at 

vector  of  natural  frequencies  (Hz) 

random  excitation  time  history 

Fourier  transform  of  f(t) 

identity  matrix 

S-i 

mass  matrix 

generalized  mass  matrix 

real  component  of  admittance,  Eq.  (6) 

imaginary  component  of  admittance,  Eq.  (6) 

frequency  of  excitation 

half -power  point  frequencies 

circular  natural  frequency 

target  mode  natural  frequency 

vector  of  modal  displacements 

vector  of  modal  velocities 

vector  of  modal  accelerations 

Fourier  transform  of  {q(t)} 

matrix  of  normal  modes 


NOMENCLATURE  ( Continued) 


♦*k 

kCh  mode  shape  value  at  degree  o£  freedom  l 

4> 

estimate  of  kc^  mode  shape  vector,  elements  of 
shaker  position  coordinates 

U(t)} 

vector  of  physical  displacements 

{X(a»)> 

Fourier  transform  of  {x(t)} 

XE 

test  article  coordinate 

C 

critical  damping  ratio 

Subscripts 

k 

counter  of  modal  coordinates 

i 

counter  of  physical  coordinates 

m 

number  of  modes  in  isolation  group 

r 

number  of  physical  coordinates 

c 

target  mode 
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